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ABSTRACT 

We report the discovery of gamma-ray emission from the Galactic globular cluster Terzan 5 using 
data taken with the Fermi Gamma-ray Space Telescope, from 2008 August 8 to 2010 January 1. 
Terzan 5 is clearly detected in the 0.5-20 CeV band by Fermi at ~ 27 'a level. This makes Terzan 5 
as the second gamma-ray emitting globular cluster seen by Fermi after 47 Tuc. The energy spectrum 
of Terzan 5 is best represented by an exponential cutoff power-law model, with a photon index of 
~ 1.9 and a cutoff energy at ~ 3.8 GeV. By comparing to 47 Tuc, we suggest that the observed 
gamma-ray emission is associated with millisecond pulsars, and is either from the magnetospheres or 
inverse Compton scattering between the relativistic electrons/positrons in the pulsar winds and the 
background soft photons from the Galactic plane. Furthermore, it is suggestive that the distance to 
Terzan 5 is less than 10 kpc and > 10 GeV photons can be seen in the future. 

Subject headings: gamma rays: stars — globular clusters: individual (Terzan 5) — pulsars: general 
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1. INTRODUCTION 

Owing to the high stellar densities in globular clusters 
(GCs), GCs are effectively pulsar factories with the pro- 
duction through frequent dynamical interactions. One 
should notice that more than ~ 80% of detected millisec- 
ond pulsars (MSPs) are located in GCs. Thanks to the 
extensive radio surveys, 140 MSPs have been identified 
in 26 GCs so far0 

The radio and X-ray properties of MSPs in GCs are 
found to be rather different from those located in the 
Galactic field (Bogdanov et al. 2006; Hui, Cheng & 
Taam 2009). The difference can be possibly related to 
the complicated magnetic field structure of the MSPs in 
a cluster, which is a consequence of frequent stellar inter- 
action (Cheng & Taam 2003). Gamma-ray observations 
can lead us to a deeper insight of whether the high en- 
ergy emission processes and hence the magnetospheric 
structure are fundamentally different between these two 
populations. However, this regime has not been fully 
explored until very recently. 

With the launch of the Fermi Gamma-ray Space Tele- 
scope (Fermi), we have entered a new era of high en- 
ergy astrophysics. As the sensitivity of the Large Area 
Telescope (LAT) onboard Fermi is much higher than 
that of its predecessor, EGRET, it has already led to 
many interesting discoveries. Shortly after its opera- 
tion, the gamma-ray emission from 47 Tuc in MeV— GeV 
regime have been detected with high significance (Abdo 
et al. 2009a). Since the other types of close binaries 
in a cluster, such as catalysmic variables and low-mass 
X-ray binaries, have not yet been detected in gamma- 
ray (Abdo et al. 2009b), the gamma-rays detected from 
47 Tuc are presumably contributed by its MSP popula- 
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tion. There are 23 MSPs have been uncovered in 47 Tuc 
so far. Both X-ray and gamma-ray analysis have placed 
an upper bound of ~ 60 on its possible MSP popula- 
tion. Venter, De Jager & Clapson (2009) suggest that 
curvature radiation results directly from the MSP mag- 
netospheres should be the dominant contributor of the 
observed gamma-rays in MeV-GeV regime. Apart from 
this component, the Inverse Compton scattering between 
the surrounding soft photon field and the energetic elec- 
trons/positrons from the pulsars can also have an addi- 
tional contribution to the gamma-ray emission (Cheng 
et al. in preparation). 

While 47 Tuc is located relatively close to us, ~ 4.0 — 
4.5 kpc (Harris 1996; McLaughlin et al. 2006), Terzan 5 
holds the largest MSP population among all the MSP- 
hosting GCs. Currently, there are 33 pulsars have been 
found in Terzan 5. Based on the radio luminosity func- 
tions of the MSPs in GCs reported by Hui, Cheng & 
Taam (2010), the hidden pulsar population in Terzan 5 
is found to be much larger than that of 47 Tuc. By scaling 
the upper limit of 60 MSPs in 47 Tuc to a corresponding 
limit in Terzan 5 with the aids of the luminosity func- 
tions, we found that the number of pulsars in Terzan 5 
can be as large as ~ 200 which is ~ 3 times larger than 
that in 47 Tuc. In view of this large predicted popula- 
tion, the collective gamma-ray intensity from the MSPs 
in Terzan 5 are expected to be well above the detection 
threshold with more than one-year long LAT data. In 
this Letter, we report the search and the analysis of the 
gamma-ray emission from Terzan 5 using Fermi LAT. 



2. OBSERVATIONS AND DATA ANALYSIS 

The Fermi/hAT is a pair-production telescope that de- 
tects gamma-rays with energies between ~ 20 MeV and 
> 300 GeV (Atwood et al. 2009). It operates in an all- 
sky scanning mode which scans the whole sky in every 
3 hours. In this study, we used the LAT data taken be- 
tween 2008 August 4 and 2010 January 1. To reduce 
and analyze the data, we used the Fermi Science Tools 
v9rl5p2 package and followed the data analysis threads 
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Fig. 1. — Fermi LAT 0.5-20 GeV (left) and 10—20 GeV (right) images of a 5° X 5° region centered on Terzan 5. The white circle shows 
the 95% confidence error circle of the gamma-ray source determined by using the 0.5-20 GeV image and gtfindsrc. The white crosses are 
the optical center of Terzan 5. In addition to Terzan 5, strong Galactic plane emission is also seen. Both images are smoothed by a 0.3° 
Gaussian function. The insert figures are the associated TS maps of the 2° X 2° region. All bright Fermi sources as well as the diffuse 
Galactic and extragalactic background emission are included in the background model. The black crosses denote the optical center of 
Terzan 5. The color scale indicates the TS value corresponding to the detection significance <r ss yTS. 

provided by the Fermi Science Support Centei0. We se- 
lected data with "Diffuse" events, which have the high- 
est probability of being a photon. In addition, we fil- 
tered out events with earth zenith angles greater than 
105° to reduce the contamination from Earth albedo 
gamma-rays. The instrument response functions (IRFs), 
"P6_V3_DIFFUSE" are used. We limited our analysis 
using data between 0.5 and 20 GeV with which the point- 
spread-function is better than ~ 1 degree such that the 
contamination by the Galactic plane emission is mini- 
mized. Furthermore, the sensitivity is more uniforrrfl. 

In Figure 1, the photon map in the 0.5-20 GeV band in 
the vicinity of Terzan 5 is shown. Although strong Galac- 
tic plane emission is seen near Terzan 5, a gamma-ray 
source at the position of Terzan 5 is clearly detected as an 
isolated source at a level of ~ 27<r (see below). We used 
gtfindsrc to determine the position of the gamma-ray 
position. In order to minimize the contamination from 
nearby sources and diffuse emission, we included Galactic 
diffuse model (gll_iem_v02.fit) and isotropic background 
(isotropic_iem_v02.txt), as well as all point sources in 
the bright source catalog^ (Abdo et al. 2009b; version 
of 2010 January) within a region of interest of 15° cen- 
tered on Terzan 5. The best- fit Fermi position of Terzan 
5 is R.A.=17h48m00s, decl.=-24d48ml5s (J2000) with 
a 95% error of 0.09° (including a 40% systematic error 
based on the bright source catalog). The optical center 
of Terzan 5 is at R.A.=17h48m05s, decl.=-24d46m48s 
(J2000), which is 0.03° from the Fermi position. There- 
fore the Fermi source is consistent with the location of 
Terzan 5. In addition, it is clear from the test-statistic 
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(TS) map in Figure 1 that the source is highly significant 
in the 0.5-20 GeV band and the optical center of Terzan 
5 falls on the peak of the TS map. Also shown in Figure 1 
is the 10-20 GeV image and its TS map; the gamma-ray 
source has less than 4cr detection significance. We also 
divided the image into four parts with roughly equal ex- 
posure time. Terzan 5 is always seen with no apparent 
flaring, indicating that the gamma-rays are not from any 
kind of flares. 

We performed spectral analysis by using maximum 
likelihood technique implemented by gtlike. Like spa- 
tial analysis, we also included diffuse emission and bright 
point sources to obtain a simultaneous fit to the unbinned 
data. All bright point sources were modeled with sim- 
ple power-laws. We first fit the spectrum of Terzan 5 
with a simple power-law; the best-fitting photon index 
r is 2.5 ± 0.1 (errors are statistical only). However, as 
shown by the gamma-ray emission from 47 Tuc as well 
as other MSPs, the energy spectrum is likely to be in 
the form of an exponential cutoff power-law (Abdo et al. 
2009a, 2009c). We show the Fermi spectrum of Terzan 5 
in Figure 2 and it is clear that the energy distribution 
appears to turn over at energies above ~ 5 GeV. We then 
re-fit the spectrum with an exponential cutoff power-law 
model; the photon index is T = 1.9 ± 0.2 and the cutoff 
energy is E c = 3.8 ± 1.2 GeV. The additional cutoff com- 
ponent is statistically significant at 99.99% level. The 
resulting TS value, which indicates the significance of 
the source detection is 725 corresponding to 27c signifi- 
cance. We show the exponential cutoff power-law model 
together with the data in Figure 2. The 0.5-20 GeV pho- 
ton flux of Terzan 5 is (3.4 ± 1.1) x 10~ 8 photons cm~ 2 
s _1 , corresponding to (6.8 ± 2.0) x 10~ n ergs cm" 2 s _1 . 
Td'rBmpare with the reported 0.1-10 GeV flux of 47 Tuc 
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Fig. 2. — Fermi LAT spectrum of Terzan 5. The solid line is 
the best-fit exponential cutoff power-law model obtained by gtlike. 
The dotted line is the best-fit simple power-law model. The error 
bars arc statistical only. 

(Abdo et al. 2009a), the 0.1-10 GeV flux of Terzan 5 
is about 2 x 10~ 7 photons cm~ 2 s _1 , corresponding to 
1.2 x 10~ 10 ergs cm~ 2 s _1 . The photon flux of Terzan 5 is 
slightly below the 95% upper limit of 2.6 x 10 -7 photons 
cm- 2 s" 1 provided by EGRET (Michelson et al. 1994). 

3. DfSCUSSfON 

Using Fermi LAT data, we detect gamma-ray emission 
from Terzan 5 for the first time. The gamma-rays emit- 
ted from GCs are generally assumed to associate with 
MSPs in clusters. The origin of these gamma-rays could 
be either pulsed curvature radiation arising near the po- 
lar cap and/or in outer magnctosphcric gaps (e.g. Zhang 
& Cheng 2003; Harding, Usov & Muslimov 2005; Venter 
& De Jager 2008), or inverse Compton scattering photons 
between the relativistic electrons/positrons in the pulsar 
winds and the background soft photons (e.g. Bednarek 
& Sitarek 2007). 

Let us first assume that these gamma-rays have the 
magnetosphcric origin. It is useful to compare Terzan 5 
with 47 Tuc, which is the first globular cluster detected 
with gamma-rays (Abdo et al. 2009a). The logarithm 
of central luminosity density and core radius for Terzan 
5 and 47 Tuc are 5.06 L Q pc~ 3 and 4.81 L Q pc~ 3 , and 
0.54pc and 0.52pc respectively (Harris 1996; version of 
2003 February). However the two body encounter rate 
and metallicity of Terzan 5 are higher than those of 47 
Tuc. Since these two quantities favor the formation of 
MSPs (e.g. Ivanova 2006,2008; Hui et al. 2010), it is ex- 
pected that Terzan 5 can house more MSPs than 47 Tuc. 
If we assume that the mean spin-down power (L sc i) of 
MSPs and the conversion efficiency of gamma-ray power 
are the same for Terzan 5 and 47 Tuc, the ratio of num- 
ber of MSPs between these two clusters is simply given 
by N Ter /N Tuc ~ (L lTer /L lTuc ) ~ 25 (assuming the dis- 
tance to 47 Tuc and Terzan 5 is 4 kpc and 10 kpc, re- 
spectively). This predicted ratio is substantially higher 
than the current observed number of MSPs for Terzan 5 
and 47 Tuc, which are 33 and 23 respectively. 

It is difficult to imagine that Terzan 5 really has 25 
times more MSPs than 47 Tuc. Perhaps the mean prop- 



erties of MSPs between these clusters are not the same. 
In fact, our spectral fits indicate that both of the pho- 
ton index and cut-off energy of Terzan 5 are larger than 
those of 47 Tuc. It has been suggested that the gamma- 

1/2 

ray efficiency is L 7 cx L sd (e.g. Thompson 2005). If 
the mean spin-down power of MSPs in Terzan 5 is ac- 
tually higher than that of 47 Tuc, NTer/NTuc will be 
reduced. Currently Fermi has detected 46 gamma-ray 
pulsars and their spectral cut-off energies (E c ) have ob- 
tained (Abdo et al. 2009c). Using the published data, 

we find that E c cx L 1 // with a correlation coefficient 
about 0.51. This also suggests that the mean spin-down 
power of MSPs in Terzan 5 is larger than that of 47 
Tuc. If this is true the number ratio can reduce to 
N T er/N Tuc ~ 25(2.5GeV73.8GeU) 2 - 11. However, it 
is still inconsistent with the observed number. Another 
possible factor to reduce this ratio is the distance to the 
cluster. There are several estimates of the distance of 
Terzan 5, which are between 5.5-10.3 kpc (Cohn et al. 
2002; Ortolani et al. 2007). If Terzan 5 is located at the 
lower end of the estimate, the ratio will become ^4. It 
is therefore likely that the distance to Terzan 5 is less 
than 10 kpc. 

On the other hand, if the gamma-rays result from 
inverse Compton scattering between relativistic elec- 
trons/positrons and the background soft photons, it is 
not surprised that {L~ t Ter / L~ t Tuc) ~ 10 (assuming a dis- 
tance of ~ 6 kpc) even though Nxer /Nt U c ~ 1. It is 
because the background soft photon intensity from the 
Galactic plane at the position of Terzan 5 is roughly 10 
times of that of 47 Tuc (Strong & Moskalenko 1998). Al- 
though current Fermi data cannot differentiate these two 
models, the curvature radiation mechanism and the in- 
verse Compton scattering mechanism have very different 
predictions in higher energy range. The curvature ra- 
diation mechanism from pulsar magnetosphere can only 
produce very few photons with energy higher than 10 
GeV (e.g. Cheng, Ho and Ruderman 1986). On the other 
hand, if E c corresponds to the peak of the inverse Comp- 
ton scattering with either relic photons (or IR photons 
from the Galactic plane), then there should be two (or 
one) more peaks correspond to IR photons and star lights 
(or star lights) respectively. Therefore the spectrum can 
be extended to 100 GeV (e.g. Bednarek & Sitarek 2007; 
Cheng et al. 2010), which may be detected by MAGIC 
and/or H.E.S.S. In particular, the background soft pho- 
tons of Terzan 5 is so high that the possibility of detecting 
photons with energy > 10 GeV is very likely. 

Because of this, we attempted to search for any > 10 
GeV photons from Terzan 5. By using the LAT data 
between 10 and 20 GeV, a maximum likelihood analy- 
sis yields a TS value of 14, corresponding to a detection 
significance of 3.7o\ Visual inspection of the 10-20 GeV 
image and its TS map (see Fig. 1) also revealed a hint 
of photon excess at the position of Terzan 5. For com- 
parison, we also analyzed the Fermi LAT data of 47 Tuc 
taken in the same period; the detection significance is 
consistent with zero in the 10 — 20 GeV band. Although 
there is an indication for a possible detection in the 10-20 
GeV band, we do not claim that it is a significant detec- 
tion for Terzan 5 since the diffuse Galactic plane emission 
is strong and there are several nearby bright gamma-ray 
sources. Deeper Fermi observations in the future and im- 
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proved knowledge of Galactic diffuse background above 
10 GeV are required to verify this speculation. 

In summary, we have detected gamma-ray emis- 
sion from the GC Terzan 5 for the first time using 
Fermi/h AT. The energy spectrum can be best fit with an 
exponential cutoff power-law model. By comparing with 
47 Tuc, we suggest that the mean spin-down power of 
MSPs in Terzan 5 is higher than those in 47 Tuc if we be- 
lieve that the gamma-rays are from the magnetospheres. 
In addition, the distance to Terzan 5 may be closer (< 10 
kpc) than expected. Alternatively, the gamma-rays can 
be from the inverse Compton scattering between rela- 
tivistic electrons/positrons and the background soft pho- 
tons from the Galactic plane. To distinguish between the 



two models, one has to look into the spectrum at higher 
energies (> 10 GeV). Future deeper Fermi observations 
as well as MAGIC and H.E.S.S. observations will be able 
to resolve this issue. 
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